Cold neutron scattering study on the diffuse and phonon excitations in the relaxor 

Pb(Mgi/3Nb2/3)03 



O 

o 

(N 



(N 



C/3 



H. Hiraka^'^, S. -H. Lee^, P. M. Gehring-^, Guangyong Xu^, and G. Shirane^ 
^Department of Physics, Brookhaven National Laboratory, Upton, New York 11973 
^Institute for Materials Research, Tohoku University, Sendai 980-8577, Japan 
NIST Center for Neutron Research, National Institute of Standards and Technology, Gaithersburg, Maryland 20899 

(Dated: February 2, 2008) 



Cold neutron scattering experiments have been performed to explore the energy, temperature, 
and wave-vector dependence of the diffuse scattering and the transverse acoustic (TA) phonons in 
the relaxor Pb(Mgi/3Nb2/3)03. We have observed a weak, but definitive, diffuse scattering cross 
section above the Burns temperature Td ~ 600 K. This cross section, which is most likely caused 
by chemical short-range order, persists down to 100 K, and coexists with the much stronger diffuse 
scattering that is attributed to the polar nanoregions. A systematic study of the TA phonon around 
(1,1,0) has also been carried out. The phonon is well defined for small wave vectors q, but broadens 
markedly around q = (0.1, —0.1, 0). 
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I. INTRODUCTION 

The lead-oxide perovskite Pb(Mgi/3Nb2/3)03, or 
PMN, is one of the most interesting and well-studied 
relaxor compounds to date, primarily because of its ex- 
ceptional piezoelectric properties and enormous poten- 
tial in industrial applicationsi*^. Neutron scattering has 
played a particularly important role in the study of the 
unusual lattice dynamics and diffuse scattering observed 
in single crystals of pM]N&lS&L2ia. Because of the Q^- 
dependence of the phonon and diffuse scattering cross 
sections, most studies to date have made use of thermal 
neutrons (A ~ 2 A) in order to access Brillouin zones cor- 
responding to large momentum transfer Q. However, Xu 
et al. have performed recent experiments on PMN using 
cold neutron (A > 4 A) time-of-flight (TOF) techniques 
that have succeeded in observing the diffuse scattering 
in the low-Q (1,0,0) Brillouin zonei&. Other cold neu- 
tron studies have also reported diffuse and low-energy TA 
phonon measurements made at low-Q ^^■-'^^j— . 

In this paper we report cold neutron scattering experi- 
ments on PMN using triple-axis spectroscopic techniques. 
Our goal was to explore the static and dynamic relaxor 
properties of PMN through measurements of the diffuse 
scattering and the transverse-acoustic (TA) phonons at 
small reduced wave-vector q in the (1,0,0) and (1,1,0) 
zones. The diffuse scattering in PMN is believed to be di- 
rectly associated with the formation of polar nanoregions 
(PNR) at the Burns temperature (T^ - 600 K^*), and it 
increases rapidly with coolin^^S. This diffuse scattering 
keeps increasing through the ferroelectric transition tem- 
perature {Tc — 210 KSiiSii&). Careful examination con- 
firms the presence of a very weak, but definitive, diffuse 
scattering cross section above Td- The high-temperature 
contours of this scattering are quite different in shape 
from those observed at low temperaturesSiSiiSiiLiSiiS. 
These contours persist to lower temperatures, remaining 
unchanged even after the much stronger low-temperature 
diffuse scattering sets in. The high-temperature contours 



are therefore very likely related to the atomic shifts pro- 
duced by the underlying chemical short-range order in- 
trinsic to PMN. 

The use of cold-neutron wavelengths limits our mea- 
surements to the (1,0,0) and (1,1,0) zones. Moreover, 
the TA phonon measurements are only possible in the 
(1,1,0) zone because the dynamical structure factor of 
the TA phonon is nearly zero at (1, 0, 0) (see Tabled). A 
well-defined low-energy TA phonon mode is observed at 
the smallest q, or ^ = 0.035, near (1,1,0) over a wide 
range of temperatures, where q — {(, — C,0). During the 
course of our experiments, Stock et al. reported a dra- 
matic line broadening of TA phonon for Tc < T < Td in 
the same (1,1,0) zone at C = 0.1 or highecSS.. Therefore, 
we have extended our cold-neutron measurements up to 
C = 0.1 for temperatures around Td as well as below Tc- 
Our results overlap and agree with their thermal neutron 
data. 



II. EXPERIMENTAL DETAILS 

Our experiments were carried out on the cold neutron 
triple-axis spectrometer SPINS located at the NIST Cen- 
ter for Neutron Research. The PMN sample, grown at Si- 



TABLE 1: Calculated structure factors for PMN. For \Fsoft'^ 
the parameter S2/S1 is fixed to be 1.5 according to Ref. 
See the text for 
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FIG. 1: Temperature dependence of the elastic diffuse scat- 
tering intensity measured at = 0.05 and 0.1. Data are 
normalized at 100 K. The solid and broken lines are guides to 
the eyes. Open diamonds show the results from TOF mea- 
surements—. 



mon Frascr University in Canada, has a mass of 4.8 grams 
and is the same high-quality single crystal that was stud- 
ied in previous reports^. The sample has a room- 
temperature lattice parameter of 4.04 A, so that 1 rlu 
(reciprocal lattice unit) = 1.56 A~^. The crystal [001] 
axis was mounted vertically so that data were collected in 
the (ft., k, 0) scattering plane. The initial (final) neutron 
energy Ei (Ef) was fixed at 4.5 meV (A = 4.26 A and 
k ~ 1.47 A^^) during elastic (inelastic) neutron mea- 
surements, and a Be filter was placed before (after) the 
sample to remove higher order neutrons. The (0, 0, 2) re- 
flection of highly-oriented pyrolytic graphite crystals was 
used to monochromate and analyze the incident and scat- 
tered neutron energies, respectively. The beam coUima- 
tions employed were guide-80'-80'-open. The measured 
instrumental energy resolution for this configuration is 
2/Ves = 0.24 meV as shown in the inset of Fig. and 
the background level is about 4 counts per minute. 

III. ELASTIC DIFFUSE SCATTERING 

Motivated by recent TOF measurements on PMNi^i, 
we have performed a systematic investigation of the dif- 
fuse scattering around both of the (1,0,0) and (1,1,0) 
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FIG. 2: Thermal evolution of the elastic diffuse scattering 
around (1,1,0) along the intensity "ridge" direction (thick 
line, inset). 



Bragg peaks. The temperature dependence of the dif- 
fuse scattering intensity is plotted in Fig. ^ at several 
different Q = G -I- q positions, where G represents ei- 
ther the (1, 1,0) or (1,0,0) reciprocal lattice vector and 
q = (C, — C, 0). All intensities have been normalized to 
100 at 100 K. The open-symbol data were measured at 
the same C = 0.05. We first note that the thermal evo- 
lution of the diffuse scattering intensity is very similar 
around different Bragg peaks. Another interesting point 
is that while the diffuse scattering intensity is visible just 
below Td, it grows rapidly below ~ 400 K. This result 
is in good agreement with previous high-resolution TOF 
measurements, which are shown by the open diamonds 
in Fig. E*' • In the TOF study the instrumental energy 
resolution {Ei = 2.7 meV and 2/Ves — 0.085 meV) was 
about three times better than that of the SPINS experi- 
ment. The consistency between these independent mea- 
surements is significant because it shows that the data do 
not change with improving energy resolution. Thus our 
data render an accurate portrayal of the intrinsic tem- 
perature dependence of the elastic diffuse scattering. At 
larger q, such as ^ = 0.1, which is shown by the solid cir- 
cles in Fig. n the diffuse scattering intensity takes off at 
a slightly higher temperature, so the whole curve (broken 
line) is shifted toward higher temperatures. However the 
overall behavior is still the same, i.e. the diffuse scatter- 
ing intensity is always first detected at or just below Td, 



and increases monotonically with cooling. 

Typical q-profiles of the elastic diffuse intensity mea- 
sured around the (1,1,0) peak are shown in a semi-log 
plot in Fig. 121 The intensity profiles are well described by 
a Gaussian function (Bragg peak) plus a Lorentzian func- 
tion (diffuse scattering) and a flat background. The fits 
are shown as the solid lines in Fig. |21 With increasing 
temperature, the intensity of the diffuse scattering de- 
creases, and the width of the Lorentzian increases, indi- 
cating a decrease in the correlation length. This temper- 
ature variation agrees with that measured in the (1,0,0) 
zone along the same [110] directioniSi (broken line, in the 
inset of Fig. 12). Note that the diffuse intensity becomes 
very weak even at 500 K, which is still 100 K below T^. 

Intensity contours of the diffuse scattering measured 
at 300 K and 650 K are shown in Fig. |3 At 300 K, 
which is above Tp, but below T^, the scattering is strong 
in both zones, and in agreement with previous studies 
^iS. The scattering intensity forms a "butterfly" pat- 
tern around the (1,0,0) Bragg peak, with ridges extending 
along the [1 ±10] directions. The scattering at (1, 1, 0) 
is ellipsoidal in shape, and extends along the [1 1 0] di- 
rection. These results are consistent with previous x- 
j.g^j,i7ji8ji9 g^jjj neutronSiSiiS measurements. On the other 
hand, a weak pattern of diffuse scattering is still visible 
around both of these Bragg peaks at 650 K, just above 
Td [Figs. O (c) and (d)]. This high-temperature scatter- 
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FIG. 3: Intensity contours of the elastic diffuse scattering 
well below Td (left column) and just above Td (right column) , 
around the (1, 0, 0) and (1,1,0) Bragg peaks. 
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FIG. 4: Temperature variation of diffuse scattering contours 
near (1,1,0). The background level is ~ 20 counts per 10 sec. 
Tuned high-temperature contours measured at 650 K in (d) 
(1,1,0) and (e) (1,0,0) zones. 



ing exhibits a geometry that is quite different from that 
observed at 300 K. The intensity of the high-temperature 
diffuse scattering at (0.9, 0.9, 0) in Fig.|31(c), for example, 
is less than 0.1 % of that of the low-temperature diffuse 
peak at (1, 1, 0). 

The presence of a diffuse cross section above Td is very 
surprising. Therefore, we did an extra background check 
using a single crystal of SrTiOa (~ 2 cc) and the identical 
configuration as used with PMN. No meaningful signal 
was observed near (1, 0, 0) and (1, 1, 0) within the exper- 
imental accuracy. 

In Fig. 01 the intensity scale was tuned to better il- 
lustrated the contrast of the high-temperature diffuse 
contour (HTC). The intensity level of the HTC is much 
weaker than that of the low-temperature diffuse peak as- 
sociated with the PNR formation as listed in Table ITTl 
However, because of the clean window and the efficiency 
of the spectrometer, the signal to noise ratio reaches 
about 4 even in the weak HTC. The HTC have distinctly 
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TABLE 11: Observed peak intensity (counts/min). 
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different shapes than the low temperature diffuse inten- 
sities in both (1,1,0) and (1,0,0) zones. They seem 
to peak at incommensurate positions around (1.1,0,0) 
and (0.9,0,0) in the (1,0,0) zone, and (1.1,1.1,0) and 
(0.9,0.9,0) in the (1,1,0) zone, while the low tempera- 
ture diffuse scattering intensities always reach the maxi- 
mum at the zone center q = 0. The intense spots at the 
center of the contour plots are scattering intensities from 
the Bragg peak. On cooling, around 450 K, the HTC 
can be clearly seen, superimposed with the increasingly 
strong low temperature diffuse cross section [see Figs. 0] 
(a) - (c)]. The weak HTC are further seen down to even 
lower temperatures in Figs. 01 (a) and (b) , and virtually 
unchanged below 650 K. Details on this high tempera- 
ture diffuse cross section will be further discussed in the 
discussion section. 



IV. PHONON AND BROAD DIFFUSE 
SCATTERING 

To date the transverse acoustic (TA) phonons in PMN 
have mainly been studied using thermal neutrons. All 
measurements are in good agreement in the (2, 0, 0) and 
(2, 2, 0) zonesi^i^ where diffuse scattering is very weak 
(see Table U). Naberezhnov et al^ discovered that the 
TA-phonon starts to broaden at Td- Later on, Wakimoto 
et al.^ foimd this slight broadening disappears at Tq. 
However, pronounced TA-phonon anomalies were ob- 
served in (3, 0, 0) and (1, 2, 2) zonea^iSii, where the diffuse 
cross section is large (see Table QJ. These anomalies were 
interpreted as the results of different types of phonon 
interactions. Our current measurements are limited in 
small Q zones as a result of the scattering geometry for 
cold neutrons, and are centered around (1, 1, 0). (1, 1, 0) 
is a good zone to study both TA phonon and diffuse scat- 
tering, because interaction from the soft transverse-optic 
(TO) phonon must be very weak due its small cross sec- 
tion. 

The fine energy resolution intrinsic to cold neutron ex- 
periments {Ef ~ 4 meV and 2/^6^ ~ 0.2 meV for typ- 
ical triple-axis measurements) has the great advantage 
of enabling the measurement of low-energy excitations 
at small q. To clarify the relationship between the dif- 
fuse scattering and the TA phonon in PMN, the dynamic 
scattering function S{q,uj) was studied along the strong 
diffuse intensity ridge q = (C, 0) for 0.035 < C < 0.1. 



The TA-phonon structure factor at (1,1,0) is much larger 
than that at (1,0,0), as demonstrated in Tabled The 
diffuse scattering structure factor listed there 

is calculated from the atomic shifts derived from neutron 
diffuse scattering measurements by Vakhrushev et al.^ 

Energy spectra measured at C = 0.035 for several dif- 
ferent temperatures are shown in Fig. |S1 A well-defined 
TA phonon mode is observed at the smallest q over a 
wide range of temperatures in spite of the relatively small 
value of Q and the small resolution volume. At 400 K 
the energy spectrum shows an interesting and definite 
overlap of the TA phonon with the central diffuse scat- 
tering intensity. This overlap could be the result of a 
softening of the TA phonon energy and the prominent 
tail of the elastic central peak. This suggests that ex- 
tra spectral weight exists in the central peak in addition 
to the resolution-limited elastic diffuse intensity. A small 
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FIG. 5: Temperature variation of constant-Q scans measured 
at small q near (1,1,0). Fits are described in the text. 
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FIG. 6: q variation of constant-Q scans measured at 300 K 
near (1,1,0). 
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FIG. 7: Constant-Q scans at = 0.1 for three temperatures. 
The incoherent scattering is shown by a thin line only in (a) . 



portion of the diffuse tail still remains at 100 K, but com- 
pletely disappears at 600 K. From these data we conclude 
that no elastic diffuse scattering is present at 600 K at 
this q. 

Energy spectra measured at 300 K and for q up to 
C, — 0.1 are shown in Fig. The overlap of the TA 
phonon with the broad central peak is still evident in 
this region. It is interesting to note that the TA phonon 
broadens markedly with increasing q at 300 K, as is quite 
apparent at C = 0.1. In order to overlap with the thermal 
neutron measurements of Stock et alm^, we have carried 
out further measurements at C = 0.1 for 600 K (~ Td) as 
well as 100 K Tc). Figure [7| traces the evolution of 
the TA phonon at (1.1, 0.9, 0). The broadening reaches a 
climax at the intermediate temperature 300 K. At T = 
100 K, which is well below Tc, the TA phonon at C = 0.1 
almost recovers its original line shape. 



In the course of the current data analyses for energy 
spectra typically shown in Figs.|5l-|7| we used the follow- 
ing scattering function to fit : 



S{q,uj) = SriD + Si 
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Here, Srw represents the resolution-limited elastic dif- 
fuse peak and it is expressed as a Gaussian function with 
the full width of 2rres as previously done^iiSii^. SbrD 
is the broad diffuse scattering prominent for intermedi- 
ate temperatures Tc < T < Td- St a describes the TA 
phonon scattering. We model SbrD and Sta in the fol- 
lowing way: 
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FIG. 8: The energy-integrated intensity of the diffuse scat- 
tering as a function of temperature. The shaded region rep- 
resents the contribution from the broad central component. 
The resolution-hmited elastic component is shown by the solid 
line. All curves are guides to the eye. 



In this model we regard the broad diffuse scattering as a 
part of intrinsic cross section for PNR. 

These functional forms in Eqs. (Q) - ® fit the data 
quite well as shown by the solid lines in Figs.|Sl-[3 where 
the phonon and the broad central diffuse component are 
depicted by the chained and dotted lines, respectively. 
Fitting parameters at C = 0.035 and 0.1 for three tem- 
peratures are listed in Table IITTI \Fta\ is derived from 
the scale factor for Eq. |21 after removing the prefactor 
of phonon intensity (T/uj^j^). The peculiar broadening 
of the TA phonon is supported numerically in Fta at 
C = 0.1 for Tc < T < Trf. In contrast, the TA-phonon 
mode at C = 0.035 is well behaved for a wide range of 
temperatures with the nearly temperature-independent 
linewidth and the dynamical structure factor. 

The temperature dependence of the energy-integrated 
intensity of the diffuse scattering at C = 0.035 {Idiff m 
Table IIII|I is shown in Fig. |S1 The broad central com- 
ponent measured at the smallest q peaks around 400 K, 



TABLE III: Representative fitting parameters without res- 
olution corrections. Idiff and \Fta\ are related to the scale 
factors in Eqs. Q and (O, respectively. 
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reaching about 10 % of the low-temperature resolution- 
limited elastic value (shaded region in Fig. (S)). Below 
400 K, which is far below Tj,, the sharp elastic compo- 
nent grows quickly and becomes dominant [see Fig.|Sl(b)] . 
As reported in previous papers^, the diffuse scattering 
intensity, corresponding to the formation of PNR, begins 
to develop around Td when observed with thermal neu- 
trons as schematically shown in Fig.|Sl(b). The difference 
between the thermal and cold neutron results can then 
be naturally attributed to differences in the experimental 
resolution functions. 



V. DISCUSSION 

The most important new result of the current exper- 
iment is the observation of the HTC in Figs. 0] (d) and 
(e). The scattering pattern is longitudinal to Q. Also 
the intensity at larger Q is stronger than that at smaller 
Q [Fig- El (e)], consistent with diffuse scattering result- 
ing from atomic shifts, where / cx |Q • e|^. Theoreti- 
cal studiesi^i^ have suggested that short-range chemical 
order between Mg^+ and Nb^+ ions at the perovskite 
B-sitcs should generate local electric field that leads to 
shifts of the A-site Pb^+ ions. This short-range chemi- 
cal order is believed to be temperature insensitive up to 
T ~ 1000 K. We speculate that this broad incommen- 
surate cross section in the HTC reflects atomic shifts 
that modulate in the longitudinal direction, caused by 
the short-range chemical order. These atomic shifts are 
also short-range correlated, resulting in the diffuse type 
scattering pattern around the incommensurate positions. 
Finally, the shape of the HTC looks somewhat conjugate 
with that of the low-temperature strong diffuse cross sec- 
tion. This may be more than merely a coincidence. It is 
likely that the diffuse contours from the PNR grow out 
of the HTC with cooling, but this is just a conjecture at 
present. 

In the diffuse-phonon-profile analysis of the (1,1,0)- 
zone data, the additional broad diffuse cross section of 
SbrD was interpreted as the part of intrinsic dynamics of 
PNR. It is the same approach employed by Gvasaliya et 
aLi^ii^ This broad Lorentzian cross section may contain 
intensity transferred from the TA phonon by coupling 
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to the resolution-limited elastic peak^i. So far we could 
not obtain a unique set of parameters by such the mode 
coupling. The cold neutron results reported here and by 
others appear to provide a complete picture of the 
low-q scattering near (1, 1, 0) and (1, 0, 0). However, pre- 
vious thermal-neutron measurements around (3, 0, 0) and 
(1,2,2) reported complex cross sections with conflicting 
interpretations ^i&L^. Thus, further neutron measure- 
ments, with high resolution in both (1,0,0) and (3,0,0) 
zones, are vital for the understanding of the PNR dy- 
namics in PMN. 
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